First-Principles Study on Electron Conduction in 
Sodium Nanowire 



Yoshiyuki Egamif, Takashi Sasakij, Tomoya Ono| and Kikuji 
Hirosef 

t Department of Precision Science and Technology, Osaka University, Suita, Osaka 
565-0871, Japan. 

X Research Center for Uhra-Precision Science and Technology, Osaka University, 
Suita, Osaka 565-0871, Japan. 

Abstract. We present detailed first-principles calculations of the electron-conduction 
properties of a three-sodium-atom nanowire suspended between semi-infinite 
crystalline Na(OOl) electrodes during its elongation. Our investigations reveal that 
the conductance is ^ 1 Gq before the nanowire breaks and only one channel with the 
characteristic of the 3s orbital of the center atom in the nanowire contributes to the 
electron conduction. Moreover, the channel fully opens around the Fermi level, and 
the behavior of the channel-current density is insensitive to the structural deformation 
of the nanowire. These results verify that the conductance trace as a function of the 
electrode spacing exhibits a flat plateau at ~ 1 Gq during elongation. 
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1. Introduction 

Metallic nanostructures have received great attention in the last decade. In particular, 
investigations of the electron-conduction properties of nanocontacts and nanowires are 
quite important since electronic devices are becoming highly integrated and miniaturized 
as nanoelectronics and nanotechnology progress. The size and structure dependences 
of the conduction properties are presently subjects of intensive research, and so far, 
numerous experimental studies on these dependences have been carried out [1-8]. In the 
experiments, atomic-size nanowires are fabricated using scanning tunneling microscopy 
or mechanically controllable break junction techniques [9]. A quantized conductance in 
units of Go = 2e^/h is observed in the nanowires consisting of monovalent atoms, such 
as Na [1,2], Cu [1,3,5], Ag [3,4] and Au [3,5-8], where e is the electron charge and h is 
Planck's constant. 

On the theoretical side, there have been many first- principles calculations performed 
for nanowires sandwiched between electrodes [10-23], and it has been reported that 
the conductances of the Au or Na single-row nanowires are quantized in units 
of Gq. However, since in some of those studies clusters [10, 11] or structureless 
jelliums [12-17] have been substituted for the electrodes, unphysical effects from the 
electrodes are always of concern. For example, in the case of jcllium electrodes, the 
conduction properties vary with the configuration of the interface between the nanowire 
and electrodes [20-23]. In addition, although the nonequilibrium Green's function 
approach with a Gaussian basis is one of the popular techniques for investigating 
the conduction properties of nanostructures [18,19], precise descriptions of electronic 
states employing this method are difficult. For a strict study, one should employ more 
realistic models including semi-infinite crystalline electrodes. Recently, the overbridging 
boundary-matching (OBM) method [24,25], which is a novel first-principles treatment of 
electron-conduction properties through a nanostructure suspended between semi-infinite 
crystalline electrodes, was developed by Fujimoto and Hirose, and several applications 
of it have been carried out [26-29]. 

In this paper, we present a first-principles study on the electron-conduction 
properties of a three-sodium-atom nanowire suspended between semi-infinite crystalline 
Na(OOl) electrodes using the OBM method. The conductance trace as a function of 
electrode spacing is observed to exhibit a fiat plateau at ~ 1 Gq before the nanowire 
breaks. By decomposing the current density into eigenchannels [14], we also find that 
only one channel is relevant to the electron conduction. This channel arises essentially 
from the contribution of the 3s orbital of the center atom. In addition, the fiat 
conductance trace before breaking is interpreted on the basis of the results, indicating 
that the characteristic of the conduction channel is negligibly affected by the structural 
deformation of the nanowire and that its transmission remains fully open around the 
Fermi level. 

This paper is organized as follows: we describe the computational procedure in 
Sec. 2, our results are presented and discussed in Sec. 3, and we summarize our studies 
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in Sec. 4. 



2. Computational Procedure 

Our first-principles calculation method is based on the real-space finite-difference 
approach [30-32] within the framework of the density functional theory [33,34]. Figure C] 
shows the schematic view of the calculation model. The nanowire is connected to the 
square bases made of sodium atoms and all of them are suspended between the semi- 
infinite crystalline Na(OOl) electrodes. The distance between the electrode surface and 
theb well as that between the basis and the edge atom of the nanowire, is initially 

set to be ao/2, where ao (=8.1 a.u.) is the lattice constant of sodium bulk; these are 
eventually determined by structural optimization. 



* t-c. I 



Figure 1. Schematic view of calculation model. The three-sodium-atom nanowire 
is suspended between semi-infinite crystalline Na(OOf) electrodes. The rectangle 
represents the supercell employed for the structural optimization and the determination 
of the Kohn-Sham effective potential. The electrode spacing is elongated gradually. 



The calculation procedure consists of two main steps. (i) The structural 
optimizations for the atomic geometry during the elongation are implemented using 
the conventional supercell of Sqq x 3ao x (2ao -l- L^s) represented by the rectangle in 
Fig. HJ where L^s is the electrode spacing. The supercell is imposed periodic boundary 
condition in x, y and z directions. The atoms in the nanowire and bases are relaxed 
by the structural optimization while the other atoms within the electrodes are fixed. 
The nine-point finite-difference formula is employed for the derivative arising from the 
kinetic-energy operator in the Kohn-Sham equation, and the grid spacing is taken to be 
0.30 a.u. The ion cores are represented by the norm-conserving pseudopotential given 
by Troullier and Martins [35] and the exchange-correlation effects are treated by the 
local density approximation [36]. (ii) The conductance is examined using the optimized 
atomic geometries obtained in step (i). Here, the central finite-difference formula, the 
local pseudopotential, the local density approximation and a grid spacing of 0.81 a.u. 
are employed. The Kohn-Sham effective potential is determined self-consistently using 
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the conventional supercell shown by the rectangle in Fig. ^ In the calculations of the 
electron-conduction properties, the supercell is imposed the periodic boundary condition 
in X and y directions and the nonperiodic boundary condition in the z direction. We 
have clarified that the enlargement of the supercell does not significantly change the 
electron conduction property. The global wave functions for infinitely extended states 
continuing from one electrode side to the other are determined by employing the OBM 
method. The electron transmission of the system at a zero bias limit is evaluated using 
the Landauer formula [37] 

G{Ep) = GoJ2T^iEp), (1) 

i 

where G{Ep) is the conductance and Ti{Ep) is the transmission of the i-th conduction 
channel at the Fermi level Ep. The eigenchannels are investigated by diagonalizing the 
Hermitian matrix (T^T) [14], where T is the transmission matrix. 

3. Results and Discussion 

The optimized geometries of the nanowire during stretching and the density distributions 
of the electrons incident from the left-hand electrode at the Fermi level are illustrated 
in Fig. 121 The planes shown are perpendicular to the [100] direction and include 
the nanowire axis. The nanowire exhibits a bent structure below L^s = 29.04 a.u. 
(Figs. |2ta) - (c)). When elongated up to L^s = 29.72 a.u., the nanowire manifests a 
geometrical transition from a bent structure to a straight one (Figs. |2l^d) and (e)), and 
finally the nanowire fractures at L^s = 31.07 a.u. (Fig. Eff)). We depict, in Fig. El 
the conductance and the channel transmission as a function of the electrode spacing. 
The conductance is ~ 1 Go and the value does not significantly change throughout the 
entire range of the electrode spacing below Lgs = 30.39 a.u. The conductance becomes 
zero after the nanowire distortion. These results of the conductance are in agreement 
with the results of the previous theoretical and experimental studies [1,10,14]. The 
details of the conduction channels can be described in terms of the decomposition of 
the incident electron-density distribution into eigenchannels. In every case, there is only 
one channel mainly contributing to the electron conduction. 

In Figs. IHA) and (B), the electronic current- density distributions for a bent 
nanowire {L^s = 28.37 a.u.) and a straight one {L^s = 30.39 a.u.), which correspond 
to configurations (a) and (e) in Fig. 121 respectively, are shown. There is no 
marked difference between these current-density distributions; the behavior of the 
electronic current-density distributions is insensitive to the structural deformation of the 
nanowire before breaking. In order to explore the characteristic of the current-density 
distributions in detail, we examine the local density of states (LDOS) of the conduction 
channels around the center atom of the nanowire and their transmissions. In Fig. the 
LDOS and the transmission of the first channel are plotted as a function of the energy 
measured from the Fermi level for (A) the bent nanowire (Lgs = 28.37 a.u.) and (B) 
the straight one {Les = 30.39 a.u.). In both cases, the 3s orbital of the center atom 
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Figure 2. Optimized geometries of nanowire and density distributions of electrons 
incident from left-hand electrode at the Fermi level, (a) L^s = 27.69 a.u., (b) 
Les = 28.37 a.u., (c) L^s = 29.04 a.u., (d) L^s = 29.72 a.u., (e) L^s = 30.39 a.u., and 
(f) Les — 31.07 a.u. Each contour represents twice or half the density of the adjacent 
contour curves. The lowest contour represents 9.93 xlO~^" electron/bohr'^/eV. 
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Figure 3. Conductance and channel transmissions of three-sodium-atom nanowire as 
a function of electrode spacing. 



is dominant in electron conduction. This result supports the independence of current- 
density distributions from structural deformation, as shown in Fig. HJ The peaks in 
the LDOS and the transmission are interpreted by the oscillation of the transmission 
probability as found in the one- dimensional penetration problem [38]. The magnitude 
of the LDOS is mainly affected by the position of the center atom: in the case of the 
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Figure 4. Channel current-density distributions of electrons incident from left-hand 
electrode at the Fermi level for (A) bent nanowire (Les = 28.37 a.u.) and (B) straight 
nanowire (Les = 30.39 a.u.). 



bent nanowire, since the center atom is not located on the central axis of the current 
flow, the magnitude of the LDOS becomes lower compared with the case of the straight 
nanowire. In addition, the LDOS and the channel transmission are fairly constant 
around the Fermi level. Because the conductance of the nanowire is determined by 
the characteristic of the conduction channel and its transmission at the Fermi level, 
our results that the channel current-density distribution is insensitive to the structural 
deformation of the nanowire and that its transmission is nnity around the Fermi level 
are consistent with the experimental result that the conductance trace as a function of 
electrode spacing exhibits a flat plateau at ~ 1 Go during elongation. 

4. Conclusion 

We performed the first-principles calculations of electron-conduction properties of 
a three-sodium-atom nanowire suspended between semi-infinite crystalline Na(OOl) 
electrodes. The conductance trace of the nanowire exhibits a flat plateau at ~ 1 Go 
during stretching. For a more profound interpretation, we examined the current- 
density distributions of electrons incident from the left-hand electrode at the Fermi 
level and decomposed them into the atomic orbital of the center atom. There is 
no significant difference between the characteristics of the current flow of the bent 
nanowire and the straight one. Only one channel, which mainly has the character 
of the 3s orbital of the center atom, is almost fully open around the Fermi level. This 
result is consistent with the experimental conductance trace that has the flat plateau at 
~ 1 Go during elongation. In a future study, we will treat the nonlocal norm-conserving 
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Figure 5. LDOS around the center atom (upper panel) and channel transmission 
(lower panel) for the first channel as a function of the energy measured from the 
Fermi level. (A) Bent nanowire (Les = 28.37 a.u.) and (B) straight nanowire 
(Les = 30.39 a.u.). 



pseudopotential to more precisely investigate the conduction properties of the various 
nanowires. 
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